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ABSTRACT 

Thermography  is  a  viable  NDE  technique  for  the  characterization  of  corrosion  in  metallic  materials.  Thermography 
is  a  rapid,  non-contact,  wide  area  inspection  technique  that  is  easy  to  interpret  and  that  in  not  significantly  sensitive  to 
material  curvature.  We  have  developed  a  portable  system  and  have  characterized  the  sensitivity  of  the  technique.  Results 
will  be  presented  in  this  paper. 
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INTRODUCTION 

The  Navy  operates  in  the  most  corrosive  environment  of  any  of  the  DoD  services  or  commercial  aviation.  As  a 
result  corrosion  prevention,  detection  and  repair  are  of  outmost  importance  to  the  Navy.  AH  the  materials  in  our  platform  are 
designed  and  engineered  to  offer  the  maximum  corrosion  protection  without  limiting  their  structural  characteristics.  New 
paint  systems  are  being  developed  that  comply  with  EPA  regulations  while  providing  enhanced  protection  against  the 
environment.  No  matter  how  much  protection  we  provide  to  our  platforms  in  the  form  of  improved  materials,  sealants  and 
paints,  the  environment  ultimately  will  penetrate  the  coating  and  initiate  corrosion.  This  process  can  be  exacerbated  or 
mitigated  by  the  specific  mission  of  the  platform.  Low  flying  antisubmarine  aircraft  will  corrode  at  an  accelerated  rate  when 
compared  to  high-altitude  communications  aircraft. 

The  cost  of  corrosion  control  is  enormous  and  accelerating  as  our  fleet  ages.  This  is  especially  significant  in  today’s 
environment  of  reduced  budgets  and  decreasing  personnel.  With  fewer  new  acquisitions,  the  Navy  is  increasingly  being 
forced  to  extend  the  life  of  existing  platforms  beyond  their  original  design  life  to  meet  mission  requirements.  In  order  to 
maintain  fleet  readiness  and  safety,  improved  inspection  methods  are  required  that  can  detect  the  occurrence  of  corrosion 
quickly  and  reliably. 

New  techniques  are  being  developed  at  the  Naval  Air  Warfare  Center  aimed  at  the  reliable  and  rapid  inspection  of 
corrosion  in  Naval  aircraft.  One  of  such  techniques  is  pulsed  thermography.  This  is  a  wide  area  inspection  technique  which 
is  especially  suited  for  the  detection  of  hidden  corrosion.  It  requires  no  contact  media  to  perform  the  inspection  as  opposed  to 
more  conventional  techniques  such  as  ultrasonic  inspection.  This  technique  relies  on  the  thermal  gradients  that  result  from 
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Figure  1 .  (a)  Shows  a  standard  single  side  inspection  setup  with  the  IR  camera  and  the  arc  lamps  in 
the  same  side,  (b)  Shows  a  through  transmission  setup  which  in  general  provides  twice  the  depth  of 
resolution. 

the  interaction  between  the  thermal  fields  and  defects  present  in  the  structure.  This  technique  is  also  significantly  insensitive 
to  curvature  and  is  a  relatively  easy  interpretation  method. 

EXPERIMENTAL 

The  minimum  number  of  experimental  components  needed  to  perform  pulsed  thermography  are  a  heat  source  and  an 
infrared  (IR)  camera.  The  proper  choice  of  heat  source  is  paramount  for  a  successful  thermographic  experiment.  Such 
source  must  be  chosen  so  as  to  maximize  the  thermal  contrast  between  the  defect  and  the  surrounding  material.  In  some 
cases  one  can  maximize  the  thermal  contrast  by  heating  the  defect  exclusively  without  disrupting  the  surrounding  material 
such  as  when  using  microwave  energy  for  detecting  water  entrapped  inside  a  ceramic  material.  The  microwave  energy  will 
propagate  through  the  ceramic  material  without  heating  it  and  finally  will  be  absorbed  by  any  water  present  in  it.  This 
approach  is  not  always  possible  and  one  has  to  resort  to  exploiting  other  defect  properties  (such  as  geometrical  features)  to 
enhance  the  thermal  contrast.  The  most  common  energy  sources  used  to  thermally  excite  materials  are:  air  heat  guns, 
microwave  sources,  infrared  lamps  and  arc  lamps.  All  the  samples  studied  in  this  paper  were  thermally  excited  using  a  pair  of 
xenon  arc  lamps,  each  one  powered  by  a  5  KJoule  capacitor  bank  with  a  10  msec  discharge  time.  The  IR  camera  used  in  these 
experiments  was  a  Amber  Engineering  InSb  focal  plane  array  (128  x  128)  camera  with  silicon  optics  operating  in  the  3  -  5 
micron  spectral  range.  Figure  la  (single  side  inspection  setup)  shows  the  experimental  configuration  used  to  image  all 
defects.  Other  experimental  arrangements  are  possible  such  as  the  one  shown  on  Fig.  lb  (through  transmission  setup).  This 
setup  is  not  well  suited  for  field  inspections  where  only  one  side  is  typically  exposed,  but  has  twice  the  depth  of  resolution 
than  single  side  inspection  methods.  Figure  2  shows  the  camera  system  and  arc  lamps  used  to  acquire  data  for  this  paper. 

Any  material  can  be  inspected  thermograph ically  (such  as  metals,  ceramics  and  composites)  but  the  outcome  of  the 
inspection  will  depend  on  many  factors  such  as  sample  geometry,  sample  thickness,  amount  of  heat  used,  surface  emissivity, 
specific  heat,  material  density,  thermal  conductivity,  defect  depth  and  the  size  of  the  defect  and  other  parameters.  To  gain 
insight  into  the  effect  that  several  of  these  parameters  have  on  the  thermographic  process  we  fabricated  various  aluminum 
7075  T6  panels  with  embedded  defects.  One  panel  was  fabricated  with  flat  bottom  holes  while  the  others  had  water  pockets 
on  the  back.  Figure  3  shows  two  thermal  images  a  few  msec  after  exposure  to  a  heat  pulse.  The  figure  on  the  left 
corresponds  to  a  Vs  1/8”  thick  aluminum  panel  with  three  different  radii  flat  bottom  holes  (R=l/2”,  3/8”  and  1/4”  diameter). 


Figure  2:  Thermal  imaging  system  used  to  acquire  data. 


Figure  3.  (a)  Shows  a  1/8”  thick  aluminum  with  flat  bottom  holes  of  different  depth  ranging  from  25 
mil  to  100  mil  and  three  different  radii  of  R=l/2”,3/8”  and  1/4”.  (b)  Shows  a  1/32”  thick  aluminum 
panel  with  water  drops  of  different  radii  and  different  water  content. 


Flat  bottom  holes  are  a  crude  simulation  to  material  loss  that  results  from  the  corrosion  process,  but  offer  a  simple  means  of 
generating  standards  for  the  characterization  of  the  thermal  process.  The  distance  from  the  surface  of  the  panel  to  the  surface 
of  the  flat  bottom  holes  (also  referred  as  the  defect  depth  “d”)  ranged  from  25  mil  to  100  mil.  The  largest  contrast 
corresponds  to  the  defect  closest  to  the  surface  while  the  defect  furthest  from  the  surface  showed  the  smallest  thermal 
contrast.  The  thermal  contrast  is  defined  as  the  difference  of  the  temperature  above  the  defect  (Td  on  fig  3  a)  from  the 
temperature  of  a  point  away  from  the  defect  area  (such  as  Tb  on  Fig.  3).  The  center  to  center  distance  between  flat  bottom 
holes  was  set  to  at  least  2  diameters  to  minimize  inter-hole  proximity  effects. 


Water  entrapment  is  a  leading  cause  of  corrosion  and  debonding  in  honeycomb  structures.  Early  detection  of  it 
leads  to  economical  repairs.  Figure  3  b  shows  the  thermal  image  a  few  msec  after  exposure  of  one  of  three  panels  (panel 
thickness  were  t0=l/32”,  1/16”  and  1/8”)  fabricated  to  model  water  entrapment.  The  panel  shown  on  that  figure  corresponds 
to  a  t0=l/32”  thick  aluminum  panel  with  water  drops  on  the  back  of  different  diameters  (1/8”,  1/4”,  3/8”and  1/2”  diameter) 
and  with  different  amount  of  water  in  them.  Water  was  contained  on  the  back  of  the  panels  by  gluing  straws  of  different 
diameters  to  the  surface  and  by  filling  them  with  water  up  to  different  levels.  The  heights  of  water  used  were  d  =  1/16”,  1/8”, 
3/16”  and  1/4”.  It  is  clear  from  Figure  3b  that  the  amount  of  water  has  little  effect  on  the  thermal  contrast  (notice  the  small 
contrast  variation  among  drops  of  equal  radii).  The  straws  were  staggered  to  maximize  the  spacing  among  them  in  order  to 
minimize  proximity  effects.  It  is  interesting  to  note  the  opposite  nature  of  the  thermal  contrast  effects  when  imaging  flat 
bottom  holes  compared  to  when  detecting  water  entrapment.  When  imaging  flat  bottom  holes  it  is  the  thinner  region  above 
the  defect  what  ultimately  makes  this  region  appear  hotter  (white  region  on  Figure  3a).  In  contrast,  it  is  the  excess  water 
acting  as  a  heat  sink  what  makes  the  region  above  the  entrapped  water  appear  colder  (black  region  on  Figure  3b). 


RESULTS  AND  INTERPRETATION 

Figure  4  shows  the  contrast  data  of  all  the  samples  studied,  the  graph  on  the  left  corresponds  to  flat  bottom  hole 
samples  while  the  graph  on  the  right  corresponds  to  samples  with  water  pockets  on  the  back.  From  the  graph  4a  clear  that 
the  deeper  the  defect  (or  the  smaller  amount  of  mass  loss)  the  smaller  the  contrast  temperature  will  be.  Also,  as  the  diameter 
of  the  defect  gets  smaller  the  contrast  temperature  diminishes.  These  results  have  been  previously  explained1,2  using  a  simple 
lateral  heat  approximation  model.  The  main  result  of  the  model  is  that  the  contrast  temperature  of  a  sample  with  a  defect  can 
be  approximated  by 
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Where  the  parameters  are  shown  in  Figure  5  (left)  and  represent: 
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Figure  4.  (Top)  This  graph  shows  the  contrast  temperature  for  all  flat  bottom  holes  shown  in  figure  3a. 
(Bottom)  This  graph  shows  the  contrast  temperature  for  all  panels  with  simulated  entrapped  water  shown  in 
Figure  3b. 


ATcont  contrast  temperature 

Q  energy  deposited  on  the  surface  of  the  sample  per  unit  area 
pc  is  the  specific  heat  times  the  density  of  the  sample 

d  distance  from  surface  of  the  sample  to  the  defect 

t0  is  the  thickness  of  the  sample 
R  is  the  radius  of  the  defect 

This  equation  correctly  accounts  for  all  of  the  observed  experimental  results.  The  term  before  the  square  bracket 
provides  the  principal  contribution  to  the  thermal  contrast  from  which  the  following  can  be  verified: 

1 .  The  contrast  temperatures  (AT)  increases  linearly  with  the  amount  of  energy  deposited  per  unit  area  (Q). 

2.  The  higher  the  specific  heat-density  of  a  material  (pcT)  the  smaller  the  contrast  temperature  becomes  (AT  i). 

3.  The  closer  the  defect  is  to  the  surface  (d  —>  0)  the  larger  the  contrast  temperature  becomes  (AT  ->  oo). 

4.  As  the  defect  depth  approaches  the  panel  thickness  (d  -»  t0)  the  contrast  temperature  vanishes  (AT  -»  0). 

5.  For  a  given  defect  depth  d,  the  thicker  the  panel  (t0  -»  oo)  the  larger  the  contrast  temperature  (AT  t). 


The  term  in  square  brackets  represent  the  lateral  heat  flaw  effects  and  the  following  can  be  verified: 

6.  As  defects  approach  the  surface  (d  — »  0)  or  as  they  approach  the  far  end  (d  ->  t0),  the  lateral  heat  effects 

diminish  ([lateral  factor]-^  1)  and  the  contrast  temperature  approaches  AT  =  — 

pc 

7.  When  the  defect  size  decreases  (R  ->  0),  the  lateral  heat  factor  vanishes  ([lateral  factor]  0)  and  the  thermal 
contrast  disappears  (AT  0). 

8.  When  the  defect  size  increases  (R  oo),  then  lateral  isolation  is  approached  and  the  lateral  heat  effects  diminish 
([lateral  factor]  -»  1)  and  the  thermal  contrast  becomes  AT  =  —  —  -  — 

Pcld  K) 

9.  When  the  thickness  of  the  panel  becomes  very  large  (t0  oo),  the  contrast  temperature  has  the  limiting  value 
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Similar  analysis  can  be  conducted  for  the  problem  of  water  entrapment.  From  the  graph  on  Figure  4  (Top)  it  is  clear 
that  the  thicker  the  panels  are,  the  smaller  the  contrast  temperature  will  be.  Also,  as  the  diameter  of  the  defect  gets  smaller 
the  contrast  temperature  diminishes.  The  error  bars  on  the  graph  represent  the  effect  that  different  amounts  of  water  had  on 
the  thermal  contrast.  It  is  obvious  from  the  smallness  of  the  error  bars  that  the  effect  of  the  amount  of  entrapped  water  on  the 
thermal  contrast  is  minimal.  A  simple  model  was  previously  obtained3  that  used  a  simple  lateral  heat  approximation.  The 
main  result  of  the  model  is  that  the  contrast  temperature  of  a  sample  with  water  entrapped  can  be  approximated  by  the 
following  equation  if  a  contact  conductivity  is  used  instead  of  the  thermal  conductivity,  then 
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Where  the  parameters  are  shown  in  Figure  5  (right)  and  represent: 

h  contact  conductivity  between  aluminum  and  water 

k  contact  thermal  conductivity  of  aluminum 

r  specific  heat-density  ratio  (r  =  pc/pwcw) 

d  the  height  of  the  water  drop 

t0  is  the  thickness  of  the  aluminum  sample 

R  is  the  radius  of  the  water  drop 


This  equation  correctly  accounts  for  most  of  the  observed  experimental  behavior  of  pulsed  thermography  as  it  relates 
to  water  entrapment.  The  terms  before  the  square  bracket  provide  the  principal  contribution  to  the  thermal  contrast  from 
which  the  following  can  be  verified: 


Figure  5.  This  illustration  shows  the  main  parameters  used  when  modeling  the  thermographic  process 
for  mass  loss  due  to  corrosion  (left)  and  water  entrapment  (right) 


1 .  The  contrast  (AT)  increases  linearly  with  the  amount  of  energy  deposited  (Q). 

2.  The  higher  the  specific  heat-density  of  the  substrate  (pcT)  the  smaller  the  peak  contrast  (AT  i). 

3.  The  contrast  temperature  (AT)  is  a  linear  function  of  the  radius  R  for  small  values  of  the  term  in  brackets. 

4.  As  the  amount  of  water  grows  indefinitely  (d  — >  oo)  the  contrast  temperature  saturates  (AT  — >  Q/piCit02). 

CHARACTERIZATION  OF  CORROSION 

There  are  various  types  of  corrosion  present  in  metallic  structures  such  as  pitting,  exfoliation,  filiform,  galvanic  and 
crevice  corrosion.  Since  the  first  means  of  corrosion  protection  is  by  painting  the  exposed  surface  or  by  sealing  joined  parts, 
corrosion  will  always  start  as  hidden.  Some  of  the  locations  where  corrosion  can  be  found  are  under  the  paint,  around 
fasteners,  inside  lap  joints,  behind  panels  and  on  second  layers.  The  ability  to  detect  corrosion  therm ographically  will  depend 
on  the  amount  of  contrast  generated  relative  to  other  contrast  generating  non-uniformities.  Figure  6  shows  a  thermal  image 
of  an  aluminum  lapjoint  with  three  rows  of  fasteners.  The  fasteners  appear  dark  because  they  are  much  thicker  than  the 
structure  being  interrogated.  Around  several  of  the  fasteners  there  are  regions  that  clearly  appear  brighter  than  others  as 
indicated  in  the  figure.  These  regions  can  be  interpreted  as  probable  corrosion  regions,  but  due  to  the  proximity  of  the 
fasteners  and  the  fact  that  sealant  was  used,  it  is  difficult  to  determine  if  any  mass  has  been  lost.  At  the  time  of  this  paper  the 
part  could  not  be  cut  open  in  order  to  make  a  definite  assessment  of  the  defect.  In  order  to  make  an  exact  determination  of 
the  cause  for  the  hot  regions  a  better  knowledge  of  the  part  and  finite  element  calculations  are  needed. 

DETECTABILITY  REGION 

There  is  a  “rule  of  thumb”  that  is  regularly  used  to  determine  the  ultimate  depth  to  which  a  defect  of  a  given  size  can 
be  detected  by  thermal  methods.  This  rule  states  that  “those  planar  defects  located  at  a  depth  smaller  or  equal  to  their  cross 
sectional  diameter  can  be  detected  thermographically”  In  light  of  the  simple  models  presented  in  this  paper,  this  rule  of 
thumb  is  not  precise.  From  eq.  1  it  can  be  seen  that  any  planar  defect,  no  matter  what  its  location  or  cross  section,  can  be 
detected  thermographically  provided  that  enough  energy  is  deposited  on  the  sample  surfaces.  This  is  direct  consequence  of 
the  linear  dependence  of  the  thermal  contrast  on  the  amount  of  energy  deposited  on  the  surface.  There  is  of  course  a 
practical  limit  of  not  depositing  so  much  energy  so  as  to  damage  the  material. 

More  insight  can  be  obtained  into  this  matter  by  inverting  eq.  1  and  obtaining  a  relation  between  the  size  of  the  flaw 
“2R”  and  the  distance  of  the  flaw  from  the  surface  of  the  sample  “d”.  The  following  equation  is  obtained 


(3) 


where  the  function  “finv[x]”  is  the  inverse  of  the  transcendental  function  y=x1/(,‘x).  This  equation  can  be  used  to  divide  the  R- 
d  space  into  “detectable”  and  “non-detectable”  regions.  That  is,  a  defect  whose  radius  “R”  and  location  “d”  are  such  that  it 


Figure  6.  Aluminum  lapjoint  with  three  rows  of  fasteners.  There  is  clear  evidence  of  a 
hot  spots  around  some  fasteners. 
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Figure  7.  This  graph  shows  curves  that  separate  detectable  from  undetectable  regions  of 
defect  of  diameter  2R. 


falls  inside  the  curve  defined  by  equation  3  will  be  detectable  thermograph ically.  If  the  defect  radius  and  position  are  such 
that  they  fall  outside  the  curve,  then  for  those  experimental  parameters  the  defect  will  be  missed.  Figure  7  shows  boundary 
curves  for  various  aluminum  thickness.  Also  on  that  graph  is  the  “rule  of  thumb”  line.  Clearly  there  is  a  strong  discrepancy 
between  the  rule  of  thumb  and  the  previous  results. 


SUMMARY 

We  have  demonstrated  that  thermography  can  detect  mass  loss  and  water  entrapment.  We  have  presented  predictive 
formulas  that  can  assist  the  practitioner  in  determining  what  type  of  contrast  levels  will  be  obtained  for  a  particular 
experimental  arrangement.  We  have  determined  what  the  depth  of  resolution  is  for  the  detection  of  mass  loss  in  aluminum 
7075T6.  When  interrogating  real  structural  components,  these  simple  models  can  only  provide  approximate  answers.  For 
precise  determination  of  the  exact  nature  and  extent  of  the  defect  a  detailed  knowledge  of  the  part  is  required  together  with 
more  detailed  analysis  such  as  finite  element  models. 
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For  a  given  defect  depth  d,  the  thicker  the  panel  (t  — »  oo)  the  larger  the 
contrast  (AT— »  q/pcd). 
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ABSTRACT 

A  common  problem  found  in  advanced  structural  materials  is  water  entrapment.  This  problem  is  a  major  cause  of 
material  degradation.  In  metals  it  can  lead  to  corrosion  and  in  composites  it  adds  unnecessary  weight  to  the  structure  and  can 
lead  to  material  degradation  especially  after  freezing  and  thawing.  Thermography  has  been  investigated  as  a  means  of 
detecting  water  entrapment  in  metallic  structures.  A  simple  model  has  been  derived  that  accurately  describes  the  thermal 
response  of  these  structures  to  short  heat  pulses.  In  this  paper  results  on  Aluminum  panels  with  various  amounts  of  water 
entrapment  will  be  presented.  Sensitivity  relations  will  be  derived  and  validated. 

Keywords:  NDE,  thermography,  thermal  models,  calorimetric  models,  defects,  water  entrapment 

1.  INTRODUCTION 

The  detection  of  entrapped  water  behind  metallic  structures  is  a  difficult  problem  to  solve  in  a  nondestructive 
manner.  X-ray  radiography  is  commonly  used  1  to  detect  the  entrapped  water  but  the  technique  is  not  very  portable  in 
general  and  requires  special  safety  precautions.  Ultrasonic  attenuation  methods  have  been  successfully  used  but  suffer  from 
a  lack  of  lateral  resolution.  In  this  paper  experimental  thermographic  results  and  simple  thermographic  models  will  be 
presented  and  discussed. 

Simple  thermographic  models  have  been  previously  introduced  2>  3  to  successfully  model  the  surface  temperature 
evolution  of  panels  with  flat  bottom  holes.  In  those  papers  the  hole  diameter  and  depth  were  varied  to  simulate  different 
amounts  of  mass  loss  due  to  corrosion.  In  all  cases  those  simple  thermographic  models  sufficed  to  accurately  predict  the 
maximum  contrast  temperatures  observed  in  the  experiments.  This  paper  will  closely  follow  the  format  and  spirit  of  the 
previous  references.  Once  again  it  will  be  seen  that  simple  physical  arguments  go  a  long  way  in  describing  complicated 
thermal  phenomena. 

D=1/8"  D=1/4"  0=3/8"  D=1/2" 
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Fig.  I:  (Left)  Drawing  indicating  key  parameters  used  in 
study.  Straws  of  various  diameters  (ranging  from  D  =  1/8” 
to  1/2”)  were  glued  to  aluminum  panels  of  different 
thickness  “t0”.  The  straws  were  filled  with  various  amounts 
of  water  “d”.  (Right)  Shows  a  single  thermographic  frame 
some  time  after  a  heat  pulse.  The  dark  regions  in  the  image 
represent  cold  spots  due  to  the  presence  of  water  in  the  back 
of  the  panel. 


Experimental  Parameters  for  Drops  and  Panels 

Al  Thickness 

Straw  Diameter 

Water  Height 

D  =  1/8" 

d  =  0" 

D  =  1/4" 

d  =  1/16" 

D  =  3/8" 

d  =  1/8" 

D  =  1/2" 

d  =  3/16" 

d  =  1/4" 

Table  1 :  This  table  shows  the  range  of  parameters  used  in  this  study 

In  this  paper  a  simple  theoretical  model  that  contains  all  the  relevant  parameters  required  to  describe  the  thermal 
transient  process  of  aluminum  panels  with  water  entrapment  will  be  presented.  Some  of  the  parameters  used  in  the  model 
include  the  amount  of  entrapped  water,  the  diameter  of  the  region  with  entrapped  water,  lateral  heat  flow  effects,  anisotropic 
properties,  finite  thickness  effects  of  the  metal  plate  and  of  the  water  column.  In  spite  of  the  simplicity  of  the  model,  it  will 
be  seen  that  the  model  fits  the  experimental  data  remarkably  well. 

2.  EXPERIMENTAL  METHOD 

The  panels  used  in  this  study  were  imaged  using  a  commercial  pulsed  IR  NDT  system  (EchoTherm®,  Thermal 
Wave  Imaging,  Inc.).  Flash  heating  was  provided  by  2  linear  xenon  flashtubes  with  5  msec  flash  duration,  each  powered  by 
a  6.4  kJouIe  capacitor  bank.  The  system  was  equipped  with  a  256  x  256  pixel  InSb  focal  plane  array  camera  (Radiance  HS, 
Raytheon  Amber)  operating  in  the  2-5  micron  spectral  range.  Continuous  12-bit  data  was  acquired  at  a  120  Hz  frame  rate  for 
6  seconds  after  flash  heating. 

The  material  used  in  this  study  was  aluminum  7075-T6.  To  model  water  entrapment,  straws  of  different  diameter 
were  glued  to  the  back  of  the  panel  and  filled  with  water  to  various  heights.  Fig.  1  (Left)  shows  a  drawing  with  some  of  the 
key  parameters  used  in  the  model.  The  parameters  “t0”,  “D”  and  “d”  represent  the  thickness  of  the  panel,  the  diameter  of  the 
straws  and  the  height  of  water  added  to  the  straws  respectively.  Three  different  thicknesses  were  chosen  for  the  aluminum 
panels.  Four  different  straw  diameters  were  studied.  The  straws  were  filled  with  water  to  different  heights.  A  total  of  48 
(3x4x4)  different  conditions  were  studied.  Table  1  summarizes  all  the  parameters  used  in  fabricating  these  panels.  The  front 
side  of  the  panels  were  painted  with  a  flat  black  paint  to  avoid  emissivity  variation  problems.  The  actual  frame  taken  soon 
after  the  thermal  flash  was  shot  is  shown  in  Fig.  1  (Right).  The  dark  areas  in  that  photo  indicate  the  presence  of  water  in  the 
back  side  of  the  panel.  The  straws  were  staggered  so  as  to  keep  the  water  as  far  as  possible  from  one  another  minimizing 
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Fig.  2:  (Left)  This  fig.  shows  the  data  curves  require  to  generate  a  single  contrast  curve.  The  curve  labeled  “Water” 
represents  the  thermal  time  history  of  a  point  on  the  surface  of  the  aluminum  panel  directly  beneath  a  straw  with  filled 
with  water.  The  curve  labeled  “Dry”  represents  the  thermal  time  history  of  a  point  away  from  any  water.  The 
contrast  curve  was  generated  by  subtracting  the  previous  two  curves.  (Right)  This  figure  shows  three  contrast  curves 
for  three  identical  water  drops  but  for  three  different  aluminum  plate  thickness. 


Maximum  Contrast 
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Fig.  3:  (Left)  This  fig.  shows  four  contrast  curves  for  a  given  aluminum  plate  thickness  and  straw  diameter  but  for  four 
different  water  heights.  (Right)  This  fig.  shows  the  relation  between  the  “maximum  temperature  contrast”  for  three 
aluminum  plate  thicknesses  and  four  straw  diameters. 

interaction  effects  while  maximizing  the  number  straws  in  one  panel. 

3.  EXPERIMENTAL  RESULTS 

In  a  standard  experiment,  digital  data  acquisition  begins  a  few  frames  before  the  capacitor  banks  are  discharged 
through  the  Xenon  arc  lamps.  In  our  experiments  the  frame  rate  was  set  to  100  frames/sec.  A  total  of  660  frames  were 
acquired  and  the  entire  experiment  lasted  over  6.6  sec.  Fig.  2  (Left)  shows  the  entire  thermal  history  of  two  points  on  the 
surface  of  the  panel.  The  curve  labeled  “water”  was  taken  from  a  point  directly  below  and  in  the  center  of  a  straw  containing 
some  amount  of  water.  This  curve  characterizes  the  thermal  evolution  a  typical  site  that  has  some  water  over  it.  The  curve 
labeled  “dry“  was  taken  from  a  point  far  away  from  any  straw  containing  water  and  characterizes  the  thermal  history  of  a  dry 
aluminum  block.  Notice  that  the  thermal  history  curve  labeled  “water”  decays  faster  and  to  a  lower  value  than  the  thermal 
history  curve  labeled  “dry”.  This  is  due  to  the  fact  that  the  “water”  acts  like  a  thermal  sink,  thereby  producing  a  decrease  in 
the  surface  temperature.  Also  notice  that  the  “water”  curve  and  the  “dry”  curve  have  a  region  at  early  times  with  very 
elevated  effective  temperature.  This  region  of  the  plot  is  difficult  to  model  because  it  includes  two  effects.  The  first  effect  is 
the  direct  radiation  from  the  rapid  temperature  rise  in  the  surface  of  the  test  panel  due  to  the  initial  flash.  The  second  effect  is 
the  reflected  radiation  from  the  walls  of  the  shroud  into  the  camera.  This  effects  disappears  when  calculating  the  contrast 
curve. 

The  difference  of  the  “water”  curve  from  the  “dry”  curve  is  termed  the  “thermal  contrast”  curve  (shown  in  Fig.  2 
Left).  Thermal  contrast  curves  start  and  end  with  zero  temperature  since  the  initial  surface  temperature  and  final  equilibrium 
temperatures  are  uniform  throughout  the  entire  panel.  Fig.  2  (Right)  shows  three  contrast  curves  for  three  identical  water 
drops  behind  three  different  plate  thicknesses.  It  is  clear  from  that  figure  that  the  thicker  the  aluminum  plate  is  the  smaller 
the  peak  thermal  contrast.  Fig.  3  (Left)  shows  four  contrast  curves  for  four  different  quantities  of  water  behind  a  1/32”  thick 
aluminum  panel.  The  diameter  of  the  straws  was  kept  constant  an  equal  to  Vi ”  while  the  amount  of  water  was  increased.  All 
four  curves  are  very  similar  in  shape  and  the  maximum  contrast  value.  When  the  experiment  was  repeated  but  with  thicker 
aluminum  plate,  the  four  contrast  curves  were  very  similar  again.  These  results  suggest  that  there  is  some  form  of  saturation 
effect  when  the  parameters  used  for  the  experiments  lay  in  the  range  as  shown  in  Table  2.  It  is  important  to  note  that  an  ideal 
contrast  curve  is  impossible  to  achieve  since  it  requires  two  points  on  the  surface  of  the  panel  infinitely  apart.  Even  in  a 
more  practical  sense  it  is  very  difficult  to  generate  systematic  contrast  curves,  because  of  difficulties  in  selecting  an 
equivalent  reference  point  on  the  surface  of  the  panel.  Some  of  the  differences  seen  in  Fig.  3  (Left)  are  simply  the  result  of 


no  lateral 


Fig.  4:  This  Fig.  shows  a  schematic  representation  of  the  zero  lateral  flow  assumption  used  in  the  model. 


difficulties  in  finding  equivalent  points  over  all  areas  of  interest. 

Finally  Fig.  3  (Right)  shows  a  plot  of  the  maximum  thermal  contrast  as  a  function  of  the  straw  diameter  for  three 
aluminum  plate  thicknesses.  Two  general  observations  were  noted  from  this  graph.  First,  water  quantity  (or  height  in  this 
case)  had  negligible  effect  on  contrast  and  second,  thicker  aluminum  panels  resulted  in  smaller  contrast.  Each  data  point  on 
this  graph  is  the  average  to  the  peak  contrast  for  the  four  water  contents  studied  (d=l/16”,  1/8”,  3/16”  and  1/4”).  The  error 
bars  represent  two  standard  deviations  of  all  peak  contrast  measurements.  In  some  cases  the  standard  deviations  are  smaller 
than  the  spot  used  to  represent  the  data  point.  This  again  illustrates  the  fact  that  for  the  range  of  parameters  used  in  this 
study,  the  effect  of  water  height  was  insignificant.  We  will  see  later  on  how  our  simple  model  is  able  to  account  for  this 
effect.  The  larger  the  panel  thickness  the  smaller  the  contrast.  For  the  thickest  aluminum  panel  (t0=l/8”)  only  one  data  point 
could  be  obtained  due  to  the  small  value  of  the  maximum  contrast.  For  the  panel  with  thickness  of  1/16”  three  data  points 
were  obtained. 


4.  CALORIMETRIC  MODEL  (ZEROTH  ORDER  APPROXIMATION) 

A  simple  theoretical  model  (zeroth  order  approximation)  has  been  derived  and  it’s  based  on  simple  calorimetric 
arguments.  Wet  and  dry  regions  are  defined  in  this  model  and  it  is  assumed  that  no  energy  flows  between  them.  Fig.  4 
shows  a  schematic  representation  of  the  model,  where  the  two  regions  (dry  and  wet)  have  been  physically  separated  to  stress 
that  there  is  no  energy  flow  between  them.  In  the  next  section  a  more  refined  model  will  be  derived  that  takes  into  account 
lateral  heat  transfer  effects. 


By  using  simple  calorimetric  arguments  we  can  write  that  q,  =  p,(A,t0) c, Ti  and  q2  =  p^A^Jc^  +  p2(A2d)c2T2 
The  parameters  q„  A„  t0,  p,,  and  represent  the  energy  deposited  over  the  area  A,  of  an  Aluminum  panel  with  thickness  t0 
density  p,  and  the  specific  heat  c^  The  parameters  q2,  A2,  d,  p2,  and  c2  represent  the  energy  deposited  over  the  area  A2  of  a 
composite  panel  made  from  an  aluminum  panel  of  thickness  t0  and  a  layer  of  water  of  thickness  d,  water  density  p2  and  the 
water  specific  heat  c2.  It  was  assumed  in  this  model  that  the  initial  temperature  of  the  coupon  was  uniform  and  equal  to  zero 
degrees.  Finally,  if  it  is  assumed  that  the  energy  deposited  on  the  surface  of  the  sample  per  unit  area  is  constant,  i.e., 
q/A,  =  q2/A2  =  Q,  then  the  temperature  difference  (or  thermal  contrast)  between  both  blocks  Tj  -  T2  =  AT  will  be 


AT  = 


JQ _ if  P2C2’d 

P.c.  t„tp2c2-d  +  p,c, -t0, 


(1) 


This  equation  correctly  describes  some  of  the  features  of  pulsed  thermography  when  applied  to  this  problem. 

1.  The  contrast  (AT)  increases  linearly  with  the  amount  of  energy  deposited  (Q). 

2.  The  higher  the  specific  heat-density  of  the  substrate  (p1c1t)  the  smaller  the  peak  contrast  (AT  1) 

3.  As  the  water  content  decreases  (d  -»  0)  the  contrast  vanishes  (AT  -»  0). 

4.  As  the  amount  of  water  grows  indefinitely  (d  ->  oo)  the  contrast  saturates  (AT  Q/piC,t0). 

It  is  interesting  to  note  that  this  simple  formula  shows  a  strong  dependence  of  the  contrast  temperature  AT  on  the  amount 
of  water  “d”  inside  the  straw.  As  was  indicated  repeatedly  in  the  previous  section,  this  was  not  an  experimentally  observed 
fact.  It  is  anticipated  the  lateral  heat  conduction  effect  (not  included  in  this  model)  is  responsible  for  the  insensitivity  of  the 
contrast  temperature  AT  on  the  variable  “d”.  This  assertion  becomes  obvious  if  we  realize  that  the  thermal  conductivity  of 
aluminum  is  two  orders  of  magnitude  larger  than  that  of  water.  In  the  next  section  a  simple  model  that  takes  into  account 
lateral  heat  flow  effects  will  be  presented. 


5  CALORIMETRIC  MODEL  (FIRST  ORDER  APPROXIMATION) 

In  this  section  the  previous  model  will  be  modified  to  allow  for  lateral  heat  flow  effects.  A  “poor  man’s  finite  element 
approximation”  (i.e.  three  elements  only)  will  be  used.  All  the  elements  will  be  thermally  interconnected.  Fig.  5  shows  a 
schematic  representation  of  this  model.  In  this  model  it  will  be  assumed  that  the  in-plane  thermal  conductivity  and  the  out- 
of-plane  thermal  conductivity  are  different.  This  will  produce  the  most  general  results. 


In  this  model  it  will  also  be  assumed  that  all  the  energy  of  the  heat  pulse  is  absorbed  in  the  aluminum  plate  of 
thickness  “t0”.  As  a  result  of  this  heating  process,  the  temperature  “T0”  of  the  aluminum  plate  can  be  derived  from 
Q  =pl-c1*p*T0  where  it  is  assumed  that  the  initial  temperature  of  the  panel  was  zero.  The  quantity  KL  represents  the  lateral 
thermal  conductance  while  Kefr  represents  the  effective  out  of  plane  contact  thermal  conductance,  the  lateral  thermal 
conductance  can  be  expressed  in  terms  of  the  lateral  thermal  conductivity  of  aluminum  k,  (which  in  this  case,  since  aluminum 
is  isotropic,  is  just  the  thermal  conductivity  of  aluminum) 

K,  =  k,  •  A  /R 

Kefr  =  h  •  A2  (2) 


Similarly,  the  contact  conductance  can  be  expressed  in  terms  of  an  effective  contact  conductivity  “h”  between  Aluminum  and 
water.  A,  and  A2  are  shown  in  the  fig.  5  and  represent  the  surface  area  of  the  dry  region  and  the  wet  region  respectively.  A2 
can  be  written  as  A2  =  nR2  where  R  is  the  radius  of  the  straw  while  A,  will  be  assumed  to  tend  to  infinity  (A!  -»  oo).  At  is  not 
shown  explicitly  in  the  figure  but  represents  the  lateral  cross  sectional  area  and  can  be  expressed  as  At  =  27tR  t0.  The  set  of 
differential  equations  that  define  this  problem  is 


dT  A 

Pi  ‘  A(  -t0  *Cj  -  =  k,  •*— l(T2  -T,) 

dt  R 

dT  A 

Pi  '  ^2  ’  t-o  *  ci  ~~r~  -  T2)  +  h  •  A2  -  (T2  —  T2 )  (3) 

dt  R 

P2 '  A2  d  c2  •— j-2-  =  h  •  A2(T2  -T2) 


Where  T„  T2,  and  T2’  are  the  temperatures  of  the  different  blocks  as  shown  in  Fig.  5.  This  set  of  coupled  differential 
equations  can  be  written  in  matrix  form  as 
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dT  =  - 
—  =  k-RT 
dt 
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This  set  of  three  coupled  linear  differential  equations  can  be  solved  by  first  diagonalizing  the  matrix  R  .  The  second  step  is 
to  transform  the  problem  into  the  eigen-space  where  the  equations  get  de-coupled.  In  this  space  the  solutions  are  simple 

decaying  exponential  functions,  the  exponent  being  related  to  the  eigen-values  of  the  matrix  R  .  Finally  the  solution  in  the 


Fig.  5:  This  figure  shows  the  building  blocks  of  our  simple  model. 


eigen-space  has  to  be  transformed  back  into  the  original  coordinate  set.  After  performing  the  previous  steps  and  imposing 
the  following  boundary  conditions 


T(t  =  0)  = 


(6) 


where  T0  is  the  initial  temperature  of  the  aluminum  block  (which  can  be  derived  from  Q  =  p,-c,-p*T0)  the  solutions  can  be 
obtained  and  the  contrast  yields 
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where  a  =  —  — - —  and  r  =  1 — -  (don’t  confuse  the  variable  “t  =  time”  with  the  parameter  “t0  =  panel  thickness”), 

h  A2  R  p2c2  d 

The  parameters  X,  and  X2  are  ^ie  non-trivial  eigen-values  of  the  matrix  R  defined  above. 


The  maximum  or  peak  thermal  contrast  can  be  calculated  by  differentiating  Eq.  7  and  the  result  gives 
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It  can  be  shown  that  in  the  limit  when  the  lateral  heat  flow  effects  disappear,  i.e  kL  — >  0  (or  a  — >  0)  then  Eq.  8  reduces  to  Eq. 
1  as  it  would  be  expected.  In  practice  the  contact  conductivity  “h”  is  much  smaller  that  the  lateral  thermal  conductivity  of 
aluminum  “k”,  as  a  result  a  »  r  >  1 .  If  eqs.  7,  8  and  9  are  expanded  around  “a”  in  this  limit,  it  can  be  shown  that 
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It  is  important  to  stress  that  the  previous  three  equations  were  derived  in  the  limit  a  »  r,  which  happens  to  be  the  “material 
parameter  space”  region  considered  in  this  work  as  shown  in  Table  1.  As  an  example,  a  vanishing  amount  of  water  in  a 
straw,  i.e.  when  d  « t0,  would  not  fall  within  these  “material  parameter  space”  region  and  therefore  would  not  be  properly 
modeled  by  the  previous  equations.  A  desirable  property  of  eq.  1 1  is  that  it  depends  very  weakly  on  the  variable  “d”.  This 


weak  dependence  on  “d”  was  experimentally  observed  and  was  described  in  section  3.  The  variable  “d”  only  enters  into  eq. 
1 1  via  the  parameter  “r”.  If  “r”  «  “a”  then  (r/a)r/a  ->  1  and  AT  actually  becomes  an  independent  function  of  “d”. 

6.  ANALYSIS  OF  RESULTS 

Figure  6  (left)  shows  all  experimental  data  points  and  the  fit  to 
them  using  eq.  11.  It  can  be  seen  from  the  graph  that  this  simple  model 
fits  the  experimental  results  fairly  well.  The  largest  fitting 
discrepancies  happen  for  the  panel  with  thickness  1/32”.  Figure  6 
(right)  shows  the  !4”  straw  diameter  data  curves  for  three  different 
panel  thickness  and  a  fit  to  them  using  eq.  11.  It  can  be  seen  again  that 
the  fitting  results  for  the  1/32”  thick  panel  show  the  largest 
discrepancies.  The  starting  parameters  that  were  used  for  all  these  fits 
are  shown  in  table  2  and  were  obtained  from  standard  material  properties  handbooks  except  for  the  contact  conductivity. 
The  final  fitting  values  were  found  to  be  reasonably  close  to  these  starting  parameters.  The  starting  value  for  the  contact 
conductivity  was  chosen  to  be  the  thermal  conductivity  of  water  per  unit  centimeter.  The  final  fitting  value  was  found  to  be 
significantly  larger  than  this  initial  estimate. 

It  is  important  to  emphasize  some  of  the  problems  found  when  dealing  with  the  type  of  data  described  in  this  paper 
and  when  using  this  simple  model.  The  first  problem  that  one  is  faced  with  is  in  generating  the  contrast  curves.  As  was 
mentioned  in  section  three,  and  ideal  contrast  curve  should  be  generated  by  subtracting  a  reference  curve  taken  far  away 
from  any  region  of  the  panel  that  has  water  from  a  curve  taken  under  a  region  with  water.  This  requirement  was  difficult  to 
meet  because  of  the  finite  size  of  our  sample  and  because  of  the  many  regions  with  water  that  we  had  in  our  panel.  To 
overcome  this  difficulty  one  could  use  one  aluminum  panel  for  each  water  drop  geometry.  The  problem  with  this  approach 
would  be  that  a  large  number  of  panels  would  be  require  and  non-systematic  errors  would  be  introduced  to  the  data  such  as 
room  temperature  variations  and  varying  amounts  of  heat  deposited  onto  the  surface  of  the  sample.  Another  source  of 
random  errors  was  the  straw  fabrication.  Straws  were  not  perfectly  circular  and  the  diameter  as  a  result  had  random  errors. 

The  model  used  in  this  work  makes  some  drastic  approximations  such  as;  1)  a  very  small  number  of  blocks  with 
average  temperatures  to  describe  the  transient  behavior  of  the  system,  2)  temperature  independent  material  properties,  3)  no 
cooling  effects  introduced  through  radiation,  4)  uniform  initial  temperature  distribution,  5)  single  drop  on  panel  (as  opposed 
the  multitude  of  drops  actually  used  in  the  panels),  6)  transverse  conductivity  approximated  by  a  contact  conductance  and  7) 
does  not  take  into  account  the  thermal  properties  of  the  straw.  It  is  anticipated  that  many  refinements  to  the  model  can  be 
achieved  by  improving  on  some  of  the  approximations.  Despite  its  crudeness  the  model  is  robust  enough  to  fit  the 
experimental  data  accurately  enough. 


Aluminum  7075  T 6 

Water 

p,  =2.8  g/cm3 

p2  =  1  g/cmJ 

c2  =  0.84  J/(g  K) 

c2  =  4.1  J/(g  K) 

k,  =  1.29  W/(cm  K) 

h  =  8.2  mW/(cm2  K) 

Table  2:  Relevant  material  property  values 


Maximum  Contrast  and  Fit 


Fig.  6:  This  figure  shows  a  fit  (solid  lines)  using  eqs.  8  or  1 1  to  the  experimental  contrast  data  for  three 
different  plate  thicknesses.  The  error  bars  represent  the  range  in  peak  contrast  temperatures  for  the  various 
amounts  of  water  used  in  the  experiments. 


7.  CONCLUSION 


Two  simple  models  have  been  developed  (figs.  4  and  5)  that  to  a  first  order  approximation  describes  the  main 
features  of  thermal  pulse  analysis  when  applied  to  a  planar  flaw.  The  first  model,  zeroth  order  approximation  does  not  take 
into  account  lateral  heat  conduction  effects.  The  second  model  takes  into  account  lateral  heat  conduction  effects,  thickness 
effects,  water  size  effects,  density  effects,  thermal  transients.  The  model  correctly  predicts  the  relationships  between  the 
previous  parameters.  Eqs.  7,  8  and  9  are  the  main  output  of  the  mode  or  their  expansions  around  he  parameter  “a”  (eqs.  10, 
1 1  and  12)  .  This  relations  were  shown  to  model  correctly  the  time  dependence  of  the  thermal  contrast,  the  peak  thermal 
contrast  and  the  time  at  which  the  thermal  contrast  peaks. 
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